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Abstract. In this work, the out-of-equilibrium dynamics of the Kardar-Parisi-Zhang 
equation in (1+1) dimensions is studied by means of numerical simulations, focussing 
on the two-times evolution of an interface in the absence of any disordered environment. 
This work shows that even in this simple rich aging behavior develops. A 

multiplicative aging scenario for the two-times roughness of the system is observed, 
characterized by the same growth exponent as in the stationary regime. The analysis 
permits the identification of the relevant growing correlation length, accounting for the 
important scaling variables in the system. The distribution function of the two-times 
roughness is also computed and described in terms of a generalized scaling relation. 
These results give good insight into the glassy dynamics of the important case of a 
non-linear elastic line in a disordered medium. 

Keywords: Slow dynamics and aging (theory), self-afhne roughness (theory), kinetic 
growth processes (theory) 
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1. Introduction 

Due to its relevance on fundamental problems of condensed matter physics, the statics 
and dynamics of elastic manifolds are a subject of great current interest. In general, 
elastic manifolds can be related to problems such as domain wall motion in magnetic [1] 
or ferroelectric [2] materials, vortex matter in high temperature superconductors [HI HIE], 
grain boundary fluctuations in materials science [6], interface dynamics in deposition 
problems [7J, or crack propagation [8j. The main ingredients characterizing these 
problems might include the elastic energy of the internal degrees of freedom, the 
quenched disorder, and the interaction between different components of the system. 
It is known that when quenched disorder is present, these systems presents glassy 
characteristics such as disorder roughness exponents or non-stationary dynamics. In 
order to better understand such features as the slow out-of-equilibrium dynamics, it 
is important to study in the first place the properties of simpler models, i.e. systems 
without disorder and/or interactions. 

It has been shown that even the simple Edwards- Wilkinson equation, a very well 
studied model of interface dynamics, presents glassy behavior [HI [101 E]. Since this is an 
harmonic solvable model, it is possible to obtain from it a lot of information about the 
out-of-equilibrium dynamics, such as the multiplicative aging scaling of the two-times 
roughness, the inclusion of finite size equilibration, and the generalization of scaling 
distribution functions. It is important to understand how these properties, obtained for 
the simple EW equation, can be generalized to other models as the non-linear Kardar- 
Parisi-Zhang (KPZ) equation [12], different discrete models [7J, and even different 
dimensions. For these cases, however, the universal exponents characterizing static 
and dynamic properties might change, and then it becomes important to understand 
how such a change affects the glassy properties. 

In particular, the study of the non-linear contributions to the elasticity should allow 
a better understanding of the glassy properties of more realistic elastic models. In this 
sense, the study of the out-of-equilibrium dynamics of the pure KPZ equation serve as 
a first step to rationalize the aging behavior reported in [13j, where a driven non-linear 
elastic model with quenched disorder was analyzed. Moreover, it shall also give a first 
insight into how the general picture obtained for the EW case can be generalized to 
more complex, and analytically non-solvable, models. 

In addition, it is now clear how important fluctuations are to describe the out- 
of-equilibrium dynamics of complex systems [HI [T5l [T6| [T7j . The thermally induced 
distribution functions of a given quantity contain more information than the first 
moments, which resulted in the proposal of using these fluctuations to characterize 
different universality classes [151 EH]- The roughness distribution of elastic lines was 
studied in the steady state for various models on the saturation regime [TBI EH ED] , and 
also in the dynamic growth regime [21] . It was also analyzed in the context of depinning 
of elastic lines in random environments [22] and in relation to the 1/f noise [23]. Its 
generalization to non-equilibrium situations has also been recently considered |TT|. [2^]. 
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Within this context, it seems important to test the scaling of the out-of-equilibrium 
roughness distribution for the KPZ equation. 

Therefore, with the aim of analyzing the glassy properties of non-linear elastic 
models, the numerical solution of the KPZ equation in (1+1) dimensions is presented 
in this work. Two-times correlation functions, as the roughness and structure factor, 
are computed and compared with the EW results. The distribution function of the 
roughness in the aging regime is also considered. These results allow to generalize 
those previously obtained for the EW equation by properly including the KPZ universal 
exponents. The paper is organized as follows. The numerical details of the simulations, 
together with the out-of-equilibrium protocol, are presented in section [2j The glassy 
properties for the different two-times quantities are presented and discussed in section |3j 
Finally, section H] presents the conclusions. 



2. Numerical details 



In the present section the numerical details of the simulation of the out-of-equilibrium 
dynamics of the KPZ equation are presented. This equation describes, in (1+1) 
dimensions, the evolution of a one-dimensional non-linear elastic object embedded in 
a two-dimensional space in the limit of small fluctuations. The non-linear contribution 
is a first order correction to elasticity, thus allowing the characterization of the line in 
terms of the univaluated field x(z, t). The KPZ equation in (1+1) dimensions is then [12] 



dx(z,t) d 2 x(z,t) 



dt dz 



2 



+ A 



dx(z, t) 



2 



+ (l) 



dz 

where the thermal noise £ is characterized by 

(£CM)) = o, (2) 
(Z(z,t)Z(z , ,t , )) = 2T8(z-z , )8(t-t'), (3) 

with v the elastic constant, A characterizing the strength of the non-linear term, T the 
temperature of the thermal bath and (• • •) the average over the white noise £, i.e. the 
thermal average. 

The KPZ equation (JTJ) is numerically solved using the following finite difference 
approach: 

Xi (t + t ) = Xi(t) + t v [x i+1 (t) - 2xi(t) + x^t)} + + y/2ATt £(*). (4) 

Here, the internal dimension is replaced by a one-dimensional lattice of unitary space, 
which sets the length scales in the following. The noise is now a random number 
uniformly distributed in [—1/2,1/2). The particular choice of the finite difference 
representation of the non-linear term is given by 

*i = - x i(t)¥ + - x i(t)\ [*i(t) - *i-i{t)] + [Xi(t) - x^t)} 2 , (5) 

which not only gives the correct exponents in the KPZ universality class but also the 
correct prefactors and the correct fluctuations in the saturation regime [25]. This is 
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important in order to perform a direct comparison with the EW case. In this work, the 
parameters used to numerically solve the KPZ equation are L = 1024 (except for the 
computation of distribution functions where the values L = 256 and L = 64 are also used 
in order to highlight the scaling properties), t = 0.01, which sets the time units, v = 1, 
and A = 1. For the thermal average, Nt = 10 3 noise realizations were considered for the 
evaluation of two-times correlations. In order to compute the roughness distributions, 
Nt = 10 5 thermal noise realizations were used for L = 256 and L = 64, while Nt = 10 4 
were used for L = 1024. 

The out-of-equilibrium dynamics of the KPZ equation was analyzed adopting 
the usual two-times protocol used to study glassy systems [26J. Starting from a flat 
configuration, the system is first equilibrated at a given initial temperature T . This 
means that the system evolves until the stationary state is reached, corresponding 
to the size-dependent saturation regime of the roughness. Then, the temperature is 
suddenly changed to the working temperature T, and the time is set to zero. Two-times 
correlation functions are then defined in terms of the waiting time t w and the elapsed 
time At = t — t w . In the present work, the working temperature was set to T = 1 
and different initial temperatures, both smaller and larger than T, were studied. In 
particular, the values T = 0, which corresponds to a perfectly flat initial condition, and 
T = 5 were used. 

3. Results 

In the following, the results for different two-times correlation quantities are described. 
Particular attention is given to correlation functions showing a clearly non-trivial out-of- 
equilibrium generalization of the scaling properties previously found in the EW equation. 

3.1. Roughness 

The one-time roughness, which characterizes the width of the discrete fluctuating line, 
is commonly defined as [7] 



where Sxi(t) = Xi(t) — x{t) and x{t) = L^ 1 J2i=i x i(t) is t ne instantaneous center of mass 
position. In general, the time evolution of the roughness can be described through the 
Family- Vicsek scaling [27], for which the roughness can be written as w 2 (t) ~ L^f(t/t x ), 
with t x ~ L z the saturation time and the scaling function f(u) ~ u 13 for u <C 1 and 
f(u) ~ const, for u 3> 1. This defines the growth exponent /3, the roughness exponent £ 
and the dynamic exponent z, related through the scaling relation z = (//3. The values 
of these exponents are £ = 1, z — 2, and (3 — 1/2 for the EW universality class, and 
£ = 1, z — 3/2, and (3 = 2/3 for KPZ universality class, both in (1+1) dimensions. 




(6) 
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Figure 1. Two-times evolution of the roughness for (a) To = and (b) To = 5. 
The different waiting times are t w — 0, 1, 8, 64, 512, 4096, and 32768, with the arrow 
indicating how t w grows. For the two largest t w values the curves are superimposed, 
indicating that the result is stationary. The dashed line indicates the (3 = 2/3 slope. 



The two-times roughness, which is the out-of-equilibrium generalization of the one- 
time roughness ([6]), is defined as 

1 L 

w 2 (At,t w ) = -^{{5 Xl {t w + At) -dx^)} 2 ). (7) 
i=i 

This quantity measures the relative fluctuation between the line's configuration at times 
t w and t. It was recently shown that, for the EW equation with finite size L, the two- 
times roughness ages [IT] , now scaling as w 2 ~ F(t w /t x , At/t x ). The aging regime 
may be understood as the pre-asymptotic non-equilibrium regime prior to finite size 
equilibration, which takes a time t x ~ L z . When t w ^> t x , the stationary solution 
and the Family- Vicsek scaling are recovered. On the other hand, in the aging regime 
t w "C t x , the scaling function behaves as 

for At«t Xl 

F[-^^)-{ (8) 



f t 

L x b x / 



L c / a (£) for At>t, 



XI 



with (3 = 1/2 and £ = 1 the EW scaling exponents. In the case where At <C t 
the roughness scales as w 2 ~ tPf^At/ty,), with ~ c(T)W for At < t w and 

~ co(T , T)^ for At 3> t w . This scaling describes a multiplicative aging scenario[j] 
for the two-times roughness [9j [TTJ [10] , and effectively corresponds to the infinite size 
limit, L — > oo. On the other hand, when At 3> t x , i.e. in the saturation regime, the 
scaling function f2(t w /t x ) behaves as fiiu) ~ s (T, T ) for t w <C t x and /2(m) ~ s{T) 
for t w 3> t x , thus leading to a independent saturation value for the roughness which 
evolves from Sq(T, T )L^ to s(T)L^ while t^ increases. The scaling form ([8]) generalizes 

| Note that for At <c t x the roughness can be also written as w 2 ~ At f[ (At/t w ), with /{ ~ c(T) for 
At <C and /{(u) ~ cq(T ,T) for At 3> t w . Equation (|8|) emphasizes the fact that this corresponds 
indeed to multiplicative aging. 
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Figure 2. Evolution of the saturation roughness with the waiting time for Tq = 
(circles) and T = 5 (diamonds), corresponding to the numerical solution of the KPZ 
equation. The continuous dashed lines correspond to the EW equation analytical 
solution. Dotted lines indicate the t w <C t x limit L^(To + T)/12. For comparison, 
the waiting time is rescaled in terms of the saturation times, t x — L 2 /144 and 
t x = i 3 / 2 /(15/2) 3 / 2 , for the EW and KPZ cases respectively. 



the Family- Vicsek scaling to a non-equilibrium situation. It is interesting to test if this 
generalization is valid for the KPZ equation, using its respective exponents. This point 
is addressed in the following. 

In this work, the two-times roughness w 2 (At,t w ), has been numerically obtained 
from the evolution of the discrete KPZ equation, and it is shown in figure HJ Panels (a) 
and (b) correspond to initial conditions with T = < T and T = 5 > T, respectively. 
The different curves correspond to different waiting times, as indicated. After an initial 
time interval where the components of the line fluctuate independently, w 2 ~ At, the 
line becomes correlated in the longitudinal direction and starts aging. Then, ones the 
line is correlated, two different time regimes are observed, corresponding to At larger 
or smaller than t x , and both showing aging. When the waiting time is larger than the 
saturation time, aging stops and the stationary situation is reached. 

In the growth regime, for each waiting time, the curves jump from the equilibrated 
asymptote to the non-equilibrated asymptote with increasing At. Both asymptotes are 
growing as At 13 , with j3 = 2/3 the KPZ value. Therefore, a multiplicative aging scenario 
with the corresponding KPZ growth exponent is found. In the saturation regime the 
curves jump from the s (T, T )L^ asymptotic value, corresponding to t w = 0, to the value 
s(T)L^, corresponding to the stationary solution. The values of the numerically obtained 
prefactors sq and s are equal, within numerical accuracy, to those analytically found for 
the EW equation, i.e. Sq = (T + T)/12 and s = T/6. It is well known that the KPZ 
and EW continuum equations have the same steady state solution for the distribution 
function of a given configuration in the saturation regime [7]. Furthermore, it has been 
also shown that this property holds for the discretized model used here [25] . However, the 
fact that the independent saturation values obtained in this work coincide with the ones 
for the EW equation suggests that the non- stationary solution in the saturation regime 
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Figure 3. Waiting time-dependent structure factor in the saturation regime, S!^(t w ), 
for the initial conditions (a) Tb = and (b) To = 5. The different waiting times are 
t w = 0, 1, 8, 64, 512, 4096, and 32768. The dashed lines indicates the k~ 2 behavior. 
The insets show a rescaled version of the main panels, using k^S^(t w ) against k n in 
order to emphasize the crossover between the two asymptotic values. Dotted lines in 
the insets correspond to the asymptotic value To+T and 2T. In the insets, t w increases 
from right to left. 



may also be the same for both models. Figure [2] shows the t^-dependent saturation 
value for both initial conditions, and compares it with the EW equation solution 
The asymptotic limits are clearly the same, and a small difference is observed in the 
crossover region. 

These results support the generalized Family- Vicsek scaling, although the growth 
regime is too short in time to accurately test the At/t w scaling. Larger system sizes 
would be necessary for a good test using the two-times roughness. However, a better 
way to observe the growing correlation length associated with aging is to focus on the 
two-times structure factor. This quantity also reflects the scaling properties of the 
roughness, but a better scaling test is obtained, as described in the next section. 



3.2. Structure factor 

The one-time dynamical structure factor is usually defined as 

S n (t) = L(\c n (t)\ 2 }, (9) 

where c n (t) represents the Fourier modes of a given configuration at time t, and where 
the discreetness of the lattice was already taken into account with the wave vectors 
given by k n = 2irn/L, n = 1, L. This structure factor gives information on how the 
fluctuations of the modes evolve in time, containing also information on the growing 
correlation length |JJ [28] . The definition (Q should not be confused with the dynamical 
structure factor S®(t) = L(c n (t)c_ n (0)), which gives information on how a given mode 
correlates between times t — and t > 0. This dynamical structure factor presents 
non-trivial stretched exponential relaxation characteristics in the KPZ case [29J. 
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In general, the structure factor ([HD scales as S n ~ t( 1+ ^/ z g(k n /k x ), with k x ~ t~ x l z . 
The scaling function g{u) behaves as g(u) ~ const, for k n <C fc^ and g(u) ~ for 
A; n ^> A^, indicating that there exists a growing correlation length £ ~ /c" 1 ~ When 
£ ~ L, the saturation regime is reached and ~ A;« ^ 1+< '' ) is independent of the time t. 

The two-times generalization of the structure factor is defined as [TT] 

S n (At,t w ) = L(\c n (t w + At) -c n (t w )\ 2 }. (10) 

This definition allows for the characterization of the fluctuations of a given mode between 
t w and t. For the EW equation, the two-times structure factor can be written as a 
scaling function expressed in terms of the quotients At/t x and t w /t x , in analogy with 
the roughness scaling ([8]). However, the focus here will be set on how the structure 
factor behaves on the ^-dependent saturation regime, which gives clear evidence of 
the growing correlation lengths. Then, for At 3> t x , i.e. in the saturation regime, the 
structure factor is computed for different t w values as 

S™(t w )= lim L(\c n (t w + At)-c n (t w )\ 2 ). (11) 

It was recently reported for the case of the EW equation that S^(t w ) ~ 
kw^ 1 gi(k n /k w ), with k w ~ t w ^ z [llj. The scaling function gi(u) behaves as g±(u) ~ 
(T + T)u~( 1+ ^ for k n <C k w and g\{u) ~ 2Tu~( 1+ ^ for k n ^> k w . The particular wave 
vector k w precisely separates two regimes. While the regime with large wave vectors, 
k n 3> k w , is equilibrated at the working temperature, small wave vectors, k n <C k w , 
still have memory of the initial temperature. Note that in contrast with the one-time 
structure factor (Q, the two-times generalization contains information on the dynamic 
exponent in the saturation regime due to its ^-dependence. 

Figure [3] shows the structure factor S™(t w ) as a function of k n for different waiting 
times, obtained in the present work. Panels (a) and (b) correspond to T = and 
T = 5, respectively. The structure factor clearly presents the k n decay, but changes 
between two asymptotes at a given t^-dependent value. In order to emphasize this 
behavior, the insets show the same data plotted as k\ S^(t w ) against k n . The crossover 
between the two asymptotes is the kind of behavior described by the scaling function 
g\{u). Indeed, the two asymptotic values are T + T and 2T respectively, indicated 
with dotted lines in the insets. These values correspond to those previously obtained 
for the EW equation, again suggesting a strong connection between the non-stationary 
solution of the EW and KPZ equations in the saturation regime. This feature is stressed 
in figure HI where the scaling of the structure factor with k n /k w ~ k n t\!, z is shown for the 
two initial conditions. The data are also compared with the EW equation solution [TT] . 
It can be seen that although the asymptotic values coincide, the intermediate dynamic 
regime does not perfectly match. This might be related to the fact that the dynamic 
exponent is also involved in the scaling function g(u), which goes beyond both models 
having the same roughness exponent in (1+1) dimensions. 
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Figure 4. Scaling of the structure factor presented in figure[3]with t w S^(t w ) ~ 

1 / z 

g(k n tj, ). The solution of the EW equation is also shown as dashed lines. The inset 
shows the same data in the scaled form k^t w 2 ^ z S™(t w ) ~ g(k n tw")- 



3.3. Scattering function 

Another quantity which has proved to be very useful in the analysis of the non- 
equilibrium dynamics in systems of interacting particles is the incoherent scattering 
function [30j SB E2] • In the context of elastic lines it can be defined as [11] 

1 L 

C q (At,t w ) = -J2 (e iq[Sxdtw+At) - Sx ^ w)] ) . (12) 
i=i 

It has been shown that due to the gaussian character of the variables involved in the 
solution of the EW equation, the scattering function can be written in terms of the 
two-times roughness in a very simple way as 

C q (At,t w ) = e-^ w2 ^\ (13) 

which establishes a clear relation between scattering functions and diffusion-like 
correlations. In the growth regime, At <C t x , it has been already shown that the 
roughness behaves as w 2 ~ t w fi(At/t w ), with fi{u) ~ c(T)u@ for At <C t w and 
fi(u) ~ co(To,T)v,P for At ^> t w , leading to a stretched exponential relaxation of the 
scattering function, i.e. 

( e -|c(T)g 2 A^ for 

C q (At,t w ) = ^ e _ hco(ToT)q2AtP for At>> ^ (14) 

This equation involves two stretched exponentials sharing the same exponent (3 but 
different prefactors. This relation trivially holds for gaussian variables, however its 
validity for the KPZ equation must be still tested. 

Figure 0(a) presents the scattering function for the KPZ equation obtained 
numerically at different q and t w values. The scattering function presents a strong 
g-dependent saturation regime, on which the q 2 factor of the exponential competes with 
the factor coming from the saturation regime of the roughness. Therefore, at low 
q values a saturation regime can be observed, while at large q values the scattering 
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Figure 5. (a) Two-times evolution of the scattering function for T = 0. The different 
waiting times are t w = 0, 1, 8, 64, 512, 4096, and 32768 from upper to lower curves. 
For the two largest t w values the curves are superimposed, indicating that the result 
is stationary. The values of q 2 are indicated for each set of curves, (b) Comparison of 
the two-times evolution of the scattering function (black lines) with the value obtained 
from the roughness and using equation (fl3|) (circles). The waiting times used here are 
t w = 0, 8, 64, 512, and 32768, from upper to lower curves. 



function decays faster to a near-zero value. This is clearly observed for the selected 
values q 2 = 0.025 and q 2 = 0.1 in figure 0(a). The generic form of the C q curves is 
similar to the one obtained for the EW equation. In order to test the validity of ( 1131) . 
the scattering function for q = 0.05 is plotted in figure [5](b) together with the values 

e 2~ w obtained using the results for the roughness presented in figure QJa). The results 
in figure 0(b) prove that relation ([13]) also holds for the KPZ equation. One could have 
expected this relation to hold in the saturation regime, since the steady state solution 
is the same for EW and KPZ equations. However, the fact that it also holds in all 
the dynamic At and t^-dependent growth regime is absolutely non-trivial, since the 
fluctuations of the KPZ are not expected to be gaussian distributed. The same results, 
i.e. the scattering function being directly obtained from the roughness through ( TLB"]) , 
were obtained for T = 5 (not shown here). 



3.4- Roughness distribution 



It has been previously proposed that the roughness distribution function in the 
saturation regime scales with the average roughness as the only scaling parameter [20J. 
This proposition can be extended to all the dynamic regime in steady state, allowing to 
write the roughness distribution as [21] 



P{w 2 ;t) 



w 2 (t) 



w 



w 2 (t) 



t 



(15) 



Here, the fluctuating roughness w 2 should be distinguished from its average value 
w 2 (t), which explicitly contains the time dependence. The last argument in f[T5"|) states 
that scaling works for curves with the same time scale t/t x , thus accounting for finite 
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Figure 6. (a) The bare distribution function P(w 2 ) for different system sizes L, 
waiting time t w , and elapsed time At, as indicated. The initial temperature is To = 0. 
(b) Scaled roughness distributions with x = w 2 /w 2 (At, t w ). Symbols as in panel 

(a). Different sets of curves correspond to the same pair of values (t w /L z , At/L z ). The 
dashed lines correspond to the EW solution in the stationary regime [20] . 



P (to 2 ; At, t w ) = ^2j^^j * 



(16) 



size effects. This scaling relation can also be generalized to include the waiting time 
dependence by considering also the t w /t x scale [51 EI], thus leading to 

w 2 At t w 

This last scaling relation was analitically obtained for the EW equation [TT] and 
numerically obtained for a directed polymer in random media model |24j . 

Using the results obtained in the present work a test of the scaling relation (Tl6l) 
for the KPZ equation can be performed. The result is shown in figure [61 obtained using 
the roughness distribution for three system sizes, L = 64, L = 256 and L = 1024, 
and different times At and t w . In panel (a), the bare distribution functions are shown, 
while panel (b) shows the scaling function with x = w 2 /w 2 (At, t w ). These results 

correspond to T = 0. Since the difference between the selected system sizes is a factor 
four, the saturation time is scaled by a factor A z = 8. Then to keep the quotient t w /t x 
fixed, the waiting time should also be scaled by the same factor 4 2 = 8; the same holds 
for the value of At if one wants to scale with At/t x . It is shown that all the curves in 
panel (a) collapse into three sets of curves in panel (b), corresponding to different pairs 
of (t w /t x , At/t x ) values, in agreement with the scaling relation ( Tl6l) . This scaling of the 
roughness distribution function also accounts for the saturation regime, leading to the 
known stationary solution. A direct comparison with the EW results in the dynamic 
regime is not possible here because the scaling factor is different; in the present case 
it should be 4 Z = 16 with z = 2 for EW. Instead, the stationary saturation solution 
of the EW case [21]] is plotted in panel (b), with a continuous dashed line, showing 
again that the numerical solution of the KPZ equation asymptotically tends to the 
same distribution function in (1+1) dimensions [T5]. 
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4. Conclusions 

The out-of-equilibrium dynamics of the KPZ equation has been studied in detail using 
numerical simulations. It corresponds to the relaxation of a non-linear interface from a 
given initial condition and it has been shown here that a complicated glassy dynamics 
emerges, as also does for the EW equation [TT] . The system size L was incorporated in 
the analysis, thus allowing to reach finite size equilibration after a time t x . Therefore, 
the out-of-equilibrium regime described here effectively corresponds to a pre-asymptotic 
regime before equilibrium is reached. One can imagine a very large system size whose 
equilibration time t x is much larger than the observation time, making the aging regime 
the relevant observable time regime. In a sense this is expected when quenched disorder 
is taken into account, where the dynamics becomes much slower and the equilibration 
time goes beyond the observation time. Thus, the main aging characteristics described 
here are of key importance for analyzing the non-equilibrium dynamics of problems 
which include other components like disorder, external forces or line-line interactions. 

It has been shown here that the main characteristics of the out-of-equilibrium 
dynamics of the EW equation can be extended to the KPZ equation. One should only 
be careful of using the proper scaling exponents. For instance, from the scaling of the 
two-times roughness and structure factor the existence of the growing correlation length 
is revealed, which allows to rationalize the different time regimes. In terms of elapsed 
time At and waiting time t w , it is the relative value of the correlation lengths £(At) and 
£(t w ) which defines the scaling properties in the growth regime. When considering the 
saturation regime, the size of the system £(t x ) ~ L, associated to the saturation time 
t x , should also be considered. Thus, it is the competition between these three length 
scales, £(At), £(t w ) and £(t x ), and its relation to the corresponding times through the 
dynamic exponent z, which define the aging dynamics of the system. 

Results concerning the scattering function have been also presented. This 
correlation function presents a well defined stretched exponential decay, which is indeed 
commonly observed in different glassy systems [26]. Moreover, it was shown that the 
scattering function is directly related to the roughness through equation (fl3|) . This 
simple exponential relation was previously analytically obtained for the EW equation 
based on the gaussian character of the fluctuations 1 1 1 j . Although in the saturation 
regime the fluctuations in the KPZ equation are also gaussian distributed, in the growth 
regime this is not necessarily true. Thus, the fact that both scattering function and 
roughness are related in this simple way for the whole two-times dependent dynamics is 
a non-trivial result, posing new questions about the relation between scattering functions 
and displacements fluctuations out of equilibrium, a fact already pointed out in colloidal 
glass experiments [31]. 

Finally, the scaling of the roughness distribution has been studied. The scaling in 
equation (flBT) for the distribution function, which depends on the relative time scales 
At/L z and t w /L z , has been tested here and a good collapse of the data was obtained. 
This indicates once more that the correct variables which include both dynamics and 



Aging of the KPZ equation 



13 



finite size effects are the relative scales between £(At), £(t w ) and £(t x ). 

Therefore, it has been highlighted here that the scaling relations analytically found 
for the EW equation are quite robust, allowing also for a good description of the out- 
of-equilibrium dynamics of the KPZ equation. While the present work focused on the 
numerical solution of the continuum equation, it could also be interesting to test these 
ideas in discrete models of interface dynamics. Finally, in order to further test these 
scaling relations, it would also be interesting to study higher dimensions for which the 
scaling exponents are different. 
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